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ABSTRACT 

Blazar subclasses ranging from fiat-spectrum radio quasars (FSRQs) through 
low-frequency-peaked BL Lac objects (LBLs) to high-frequency-peaked BL Lac 
objects (HBLs) exhibit a sequence of increasing spectral hardness with decreas- 
ing luminosity that cannot be explained solely by orientation effects. Using an 
analytic model for the synchrotron, synchrotron-self-Compton, and Compton- 
scattered external radiation from blazar jets, we propose an evolutionary sce- 
nario that links these blazar subclasses in terms of a reduction of the black-hole 
accretion power with time. As the circumnuclear material accretes to fuel the 
central engine, less gas and dust is left to scatter accretion-disk radiation and 
produce an external Compton-scattered component in blazar spectra. This evo- 
lutionary trend produces the sequence FSRQ LBL — > HBL. Such a scenario 
may also link radio-loud AGNs with ultraluminous infrared galaxies and optical 
QSOs, if the latter constitute the high Eddington-ratio epoch of supermassive 
black hole growth, as suggested by the observed anti-correlation between radio 
and soft X-ray activity in some Galactic black-hole candidates. 

Subject headings: galaxies: active — gamma-rays: theory 
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1. Introduction 

During the past decade, significant progress has been made in understanding the physics 
of active galactic nuclei (AGN) and the connection between different subclasses of AGN. 
In particular, various lines of evidence indicate that FR II radio galaxies are the parent 
population of radio-loud quasars, with the relativistic jets of FSRQs being closely ahgned 
to our hne-of-sight (for a review, see Urry & Padovani (1995)). Orientation effects are also 
thought to unify FR I radio galaxies with BL Lac objects. The BL-Lac/FR-I subclasses have 
lower average luminosities than the FSRQ/FR-II subclasses. Among blazars, which include 
highly polarized and optically violently variable quasars, FSRQs and BL Lac objects, there 
appears to be an almost continuous sequence of properties from FSRQs through LBLs to 
HBLs. This trend is characterized by decreasing bolometric luminosities, a shift of the peak 
frequencies of their broadband spectral components towards higher values, and a decreasing 
fraction of power in 7 rays compared with lower-frequency radiation (Sambruna, Maraschi, 
& Urry 1996; Fossati et al. 1998). 

In the framework of relativistic jet models, the low-frequency (radio - optical/UV) 
emission from blazars is interpreted as synchrotron emission from nonthermal electrons in 
a relativistic jet. The high-frequency (X-ray - 7-ray) emission is thought to be produced 
via Compton upscattering of low frequency radiation by the same electrons responsible for 
the synchrotron emission. Possible sources of soft seed photons for Compton scattering 
are the synchrotron photons themselves (Maraschi, Ghisellini, & Celotti 1992; Bloom & 
Marscher 1996) or external photons, presumably dominated by the accretion disk emission 
which can enter the jet either directly (Dermer, Schlickeiser, & Mastichiadis 1992; Dermer 
& Schlickeiser 1993) or after reprocessing by circumnuclear gas and dust (Sikora, Begelman, 
& Rees 1994; Dermer, Sturner, & Schlickeiser 1997). In addition, a significant contribution 
to the soft radiation field may also be provided by infrared emission from dust in the vicinity 
of the AGN (Blazejowski et al. 2000; Arbeiter, Pohl, & Schlickeiser 2001). 

The blazar sequence has been studied by Ghisellini et al. (1998) using a large sample 
of blazar broadband spectra. They suggest that along the sequence HBL — LBL — FSRQ, 
an increasing energy density of the external radiation field leads to an increasing amount of 
Compton cooling. As a result, the maximum energy in the electron distribution decreases, 
causing the synchrotron and Compton peaks to shift to lower frequencies. This explanation 
for the blazar sequence was developed further by Georganopoulos, Kirk, & Mastichiadis 
(2001), who argue that the radiating jet plasma is outside the broad-line scattering region in 
weak sources and within it in powerful sources. Detailed model fits to the broadband spectra 
of several blazars also indicate that a decreasing contribution of the external radiation to the 
seed photon field reproduces the FSRQ — > LBL — > HBL sequence. While FSRQs generally 
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require a dominant external Compton contribution to produce the observed 7-ray emission 
(e.g., Sambruna et al. (1997); Mukherjee et al. (1999); Hartman et al. (2001)), HBLs can be 
successfully fitted with pure synchrotron self-Compton (SSC) models (e.g., Mastichiadis & 
Kirk (1997); Pian et al. (1998); Petry et al. (2000)). The LBL BL Lacertae appears to be 
intermediate between those two classes, requiring a non- negligible contribution from external 
Comptonization to reproduce its 7-ray spectrum (Madejski et al. 1999; Bottcher & Bloom 
2000) . Orientation effects certainly play a role in spectral variations among sources within a 
given subclass, or even to account for the blue blazar subclass (Georganopoulos 2000), but 
cannot explain the full range of blazar properties (Sambruna, Maraschi, & Urry 1996). 

It is uncertain whether the different blazar types are connected through an evolutionary 
sequence, or whether they constitute limited periods of AGN activity on parallel evolutionary 
paths. An evolutionary sequence has recently been proposed by D'Elia & Cavahere (2000) 
and Cavahere & D'Elia (2001) who consider the efficiency of energy extraction from the 
accretion flow onto the central supermassive black hole through magnetic fields anchored in 
an optically thick accretion disk vs. energy extraction from rotational energy of the black 
hole via the Blandford-Znajek (BZ) mechanism (Blandford & Znajek 1977). They argue 
that BL Lac objects can plausibly be powered by the BZ mechanism without a significant 
contribution from an accretion disk, while FSRQs require a dominant contribution from the 
accretion disk to power the relativistic jets responsible for the production of the nonthermal 
broadband spectra. They suggest that a gradual depletion of the central environment in 
FSRQs by accretion onto the central black hole may be responsible for a transition from the 
disk-powered jet production (in FSRQs) to the BZ-powered jet production mode (in BL Lac 
objects). Furthermore, they find that a cosmological BL Lac birth rate proportional to the 
rate of "deaths" of FSRQs (through depletion of their central environments) is consistent 
with the cosmological evolution observed for both quasars and BL Lac objects. The masses 
of the central black holes in radio quasars generally exceed ~ 3 x 10^ Mq (Wandcl 1999; Laor 
2000). Similarly detailed mass determinations have so far not been possible for the lineless 
BL Lacertae objects. The cluster environments of FR II radio galaxies and radio quasars on 
the one hand, and FR I radio galaxies and BL Lacs on the other hand appear consistent with 
an evolution of the former into the latter class of objects (Hill & Lilly 1991; Yee & EUingson 
1993; Urry & Padovani 1995). Also, the distribution of central black- hole masses in nearby 
non-active galaxies appears to extend to larger values than the distribution of radio-quasar 
black- hole masses (Laor 2000), consistent with some fraction of present-day giant elliptical 
galaxies hosting the remnants of earlier quasar activity. 

Here we expand on the suggestion of a cosmological evolutionary sequence (D'Elia & 
Cavahere 2000; Cavahere & D'Elia 2001) from FSRQs to BL Lac objects by demonstrating 
that such a sequence plausibly and self-consistently explains the trend in their spectral 
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properties. We employ a combined SSC/external Compton model for relativistic jet emission, 
as described in §2, and derive parameter values from the observed multifrequency signatures 
of FSRQs and HBLs. The evolutionary sequence connecting the respective parameter regimes 
is presented in §3. We discuss our results and suggest a possible evolutionary connection 
between radio-loud and radio-quiet AGN classes in §4. 

2. BlciZcir Model and Pcirameter Estimation 

The radiating plasma is modeled by a spherical emission region of co-moving radius Rb 
that moves outward along the axis of the jet with bulk Lorentz factor F = (1 — 
Electrons are injected into the plasma blob with a comoving power-law distribution Q{'~f) = 
Qo'y~^, for electron Lorentz factors in the range 7i < 7 < 72- The jet power Lj corre- 
sponding to the energy input of injected particles into the jet is assumed to be proportional 
to the accretion disk power Ld, as suggested by the observed proportionality between ki- 
netic jet power and narrow line emission luminosity in radio-loud AGN reported by Rawling 
& Saunders (1991) and the proportionality between core radio luminosity and broad emis- 
sion hne luminosity in radio-loud AGN found by Celotti, Padovani & GhiseUini (1997). 
In addition, Xu, Livio & Baum (1999) have found a similar proportionality between the 
5 GHz radio luminosity and the [O 111] 5007 emission line luminosities both for radio-loud 
and radio-quiet AGN. The jet power is related to the injected particle distribution through 
Lj = fj Ld = meC^ J^^ d'j 7 (5(7)- 

Electrons lose energy through adiabatic expansion, which is treated in terms of an escape 
time scale tcsc = V^b/c with rj ^ 1, and through radiative losses. These losses are generally 
dominated by synchrotron and Thomson processes, so they can be parameterized by the 
functional form 7rad = — i^o7^, where uq = AcaT (mb + Wsy + i^ext) l'irae(? ■ The terms Ms, 
lisy, and Mext represent the comoving magnetic field, synchrotron, and external radiation field 
energy densities, respectively. If Wext is dominated by reprocessed accretion-disk radiation due 
to circumnuclear gas and dust with a scattering optical depth Trepr at a characteristic distance 
it!sc from the central source, then Wext f« Trepr/ (47r -Rsc*^)- Our approach is similar to the 

one used by GhiseUini et al. (1998), but differs from their model by our explicitly accounting 
for adiabatic losses and escape. These processes lead to the formation of a low-energy cutoff 
in the electron spectrum and force us to consider two distinct radiative regimes of electron 
cooling and of the equilibrium electron spectra (see below), while the GhiseUini et al. (1998) 
approach always leads to a 7^^ electron spectrum for 1 < 7 < 71. We neglect pair processes, 
which generally turn out to be unimportant in the course of detailed spectral modeling of 
blazars. 
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Electrons injected with 7 ^ 7^ = c/lRbUo) cool rapidly prior to escape, whereas 
electrons injected with 7 -C 7cr lose only a small fraction of their energy before they escape. 
The competition between electron injection, cooling and escape establishes an equilibrium 
distribution given by 

/NO jil/lmy^ for7i<7<7^ 

Two distinct radiative regimes can be defined. If 71 < jci, most of the injected electrons 
cool inefficiently and the electron population is in the slow cooling regime. In this case, 
p = s, 7; = 71, 'jm = Icr, and n°gg = Qo^esc7mVH, where Vl is the blob volume. When 
7i > 7cr, the electron population is in the fast cooling regime, p = 2, = jcr, 1m = 1i, and 
Hg gg = Qolm^^'^^^ /[VbUo — 1)]. The synchrotron, and external-Compton radiation emitted 
by these relativistic electron populations are evaluated using the formulae given by Dermer, 
Sturner, & Schlickeiser (1997). The SSC radiation is calculated using the analytic solution 
of Tavecchio, Maraschi, & Ghisellini (1998). 

In the absence of an external Compton component, we recover the SSC model used to 
fit data from HBLs such as Mrk 501, Mrk 421, or PKS 2155-304 (e.g., Mastichiadis & Kirk 

(1997); Pian ct al. (1998); Tavecchio, Maraschi, & Ghisellini (1998); Bednarck & Protheroe 

(1999) ; Pctry ct al. (2000); Sambruna et al. (2000); Fossati et al. (2000); Kataoka et al. 

(2000) ). The parameters in an SSC model are the Dopplcr factor D = [T {1 — /5r cos ^^obs)]"^, 
the magnetic field B, Ri,, 71, 72, the injection luminosity Lj, and the injection spectral index 
s. Ill principal, the escape time scale parameter r/ could also be considered a free parameter, 
though we assume (see Kataoka et al. (2000)) that r] ^ 1, i.e., tesc ~ Rb/c. These studies 
find typical HBL parameters of Rb ~ 10^^ - 10^^ cm, 72 ~ 10^ S ~ 0.1 G, D ~ 20 - 30, and 
a rather broad range of injection luminosities Lj ~ 10^^ ergs s~^ - 10^^ erg s~^, depending 
on the various model details and the individual flares that are modeled. 

The addition of the external Compton component in FSRQs complicates parameter es- 
timation. However, since detailed modeling of blazar spectra using self-consistent, combined 
SSC -I- ERG spectral models is generally a rather time-consuming task, simple formulae 
to achieve order-of-magnitude estimates for the relevant model parameters are desirable in 
order to find a convenient starting point for more detailed modeling efforts. In the follow- 
ing, we derive such estimates, although we must caution that some of the observables, in 
particular those pertaining to the SSC component, are rather uncertain, which may intro- 
duce significant scatter in the derived analytical parameter estimates. We assume that the 
X-ray spectra of FSRQs are dominated by the SSC process, as has been argued by various 
authors on the basis of spectral fitting as well as variability considerations (e.g., Mukherjee 
et al. (1999); Hartman et al. (2001); Sikora et al. (2001)), while the 7-ray spectrum is domi- 
nated by the external Compton component. The peak photon energies and peak i/Fi, fluxes. 



-6- 



uFi, — 10~^° -F-io ergs cm~^ s~^, of the Compton and synchrotron components are determined 
from simultaneous multiwavelength observations. In the following, photon energies are given 
in dimensionless units e = hv/me(? = 10" e„, and we define = B/B^r = 2.3 x 10~^^ B{G). 
The jets in blazars are generally believed to be oriented close to the line of sight where the 
Doppler beaming factor D ^ T. Thus, for simplicity, we assume D = T in the following. 
The model parameters can then be estimated from the observables through the expressions 

cerc = — r- — — , (3) 



1 + z 

1 + Z 



esse — -r- , (4; 



(^^^)sy = T^n t "■B7cr /sp /sy , (5) 



4ndl 3 

(i^-^f^)ERC = 3 Mext7cr /sp /erc , and (6) 

L 

(^^>^)ssc = c4 Mb 7,^, (^e/sp)^ /ssc • (7) 

The subscripts "sy," "SSC," and "ERC" refer to the synchrotron, SSC, and external Compton 
components, respectively, and is the mean photon energy of the external soft photon field 
in the stationary frame of the AGN. The factor /gp in equations (5) - (7) is a normalization 
factor defined through f^djne{j)j^ = 7cr^e/sp; and /ssc in equation (7) is a correction 
factor between the uF^ peak value and the total energy output in the SSC component, 
which differ significantly due to the substantial spectral broadness of the SSC emission. 
Typically, /ssc ~ 0.1. For the more strongly peaked synchrotron and ERC components, the 
corresponding correction factor is fgy ~ /erc ~ 1- 

Combining the above estimates, and defining (II ~ 10^^ ^28 cm and it!sc = 10^^ -Rsc.is cm, 
we find 

D = 3.2^i^^^^^^=^^±^, (8) 
e*,-5 essc,-i 




^^1.4^/ -^'-^'*'-^^^ + "^^ , G (9) 

essc.-i eERC,2 



Lu r.epr = 3 X 10^^ ^^^^ (l + ^)^L,18^sV7<-5 ^-1 ^ (^Q^ 

-fsy,-10 e^i^Qo 
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i?esc = 1.2 X 10^^ 7r?5^ f 1 + " cm , (11) 




D inl7 "^28 / -^sy-10 / <^*,-5<^SSC -1 / (/sSc/O-l) , . s 

^6 = 10^^=^\ TT^^ ^5 \/ « cm, and (12) 

Note that equation (13) requires an independent estimate of 71. Assuming that the optical 
spectrum is produced through synchrotron emission from strongly cooled electrons, the in- 
jection spectral index s can be estimated through the optical spectral index ojopt from the 
relation s = 2 aopt- 

The requirement that r] = Rcsc/Rb > 1 yields a constraint on the external soft photon 
energy e* from equations (11) and (12), given by 

^ ^ 2.4 X 10-^ 6ERC.2 ( ^/K^o f /sy 

(^rT^ci28)'^' 4c,-l V^^^'-^° + ^^^C,-10; \y/{fssc/10)j 

However, we point out that the observables esse and {j^F^)ssc generally uncertain 
by at least an order of magnitude since the SSC component does generally not show up as 
a pronounced bump in the broadband spectra of FSRQs. Also the ERC peak energy eERc 
can rarely be determined to an accuracy of better than an order of magnitude. Thus, great 
caution must be used when drawing conclusions from Eq. (14). 

For the example of 3C 279, one of the best-observed EGRET-detected FSRQs, we 
can estimate typical values of the observables as egy ~ 3 x 10~^, cerc ~ 300, esse ~ 0.1, 
Fsy _io ^ 1, Ferc,-io ~ 5, and Fssc,-io ~ 0.5 (sec, e.g., Hartman et al. (2001)). For 3C 279, 
z = 0.538 and ^28 = 0.77, and we take Tj-cpr = 0.1, 71 = 100, and s = 2.2. Equation (14) 
then implies <,1.8 x 10^^ assuming rj = 1 (corresponding to ^ 1 eV or u ^ 2 x 10^^ Hz). 
Possible implications of this soft seed photon energy estimate will be discussed in §3. Other 
model parameter estimates based on these observables are D = 28, 5 = 2.2G, i?sc ~ 
1.2 X 10^^ Vi^46(7-repr/0.1) cm, Rf, = 1.7 X lO^^ ^ni, fj = 6.4 X IQ-^ LJe^ (7i/100)-°-2, implying 
Lj ~ 6.4 X 10^^ (7i/100)~°-^ ergs s~^. The high-energy cut-off of both the synchrotron and 
the Compton components are poorly known in 3C 279 and other FSRQs, so no rehable 
constraint on the maximum electron energy 72 can be derived in the FSRQ case. 

The parameter estimates for the FSRQ case quoted above are in good agreement with 
those found by GhiseUini et al. (1998) for FSRQs in general, and through more detailed 
spectral modehng of 3C 279 in particular (e.g., Hartman et al. (2001)), which indicates that 
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our analytic estimates may indeed be useful to assess a convenient starting point for more 
detailed spectral modeling of FSRQs. 

3. Evolutionary Sequence from FSRQs to HBLs 

In the previous section, we obtained typical parameter values for HBLs and FSRQs. 
These estimates indicate that the most significant differences between model parameters of 
FSRQs and HBLs are found in the values of the jet power Lj and the magnetic field B, which 

arc both found to decrease by more than an order of magnitude from FSRQs to HBLs. The 
lack of evidence for a significant contribution of an external Compton component to the 
7-ray spectra of HBLs indicates that the energy density of the external soft photon field — 
parameterized through Ld, Trcpr, and i?sc — decreases by >,3 orders of magnitude along the 
sequence of blazar types. This is in agreement with the results of Ghisellini et al. (1998). 

The decreasing energy density of the external soft photon field along the sequence FSRQ 
— > LBL HBL, in tandem with the decreasing jet power suggests that this is a consequence 
of a decreasing accretion rate along this sequence. This, in turn, is in very good agreement 
with the suggestion of D'Elia & Cavaliere (2000) and Cavaliere & D'Elia (2001) of an evo- 
lutionary sequence from FSRQs to BL Lacs characterized by such a decline in the accretion 
rate. We are thus motivated to try to reconstruct the evolutionary sequence FSRQ — > LBL 
— > HBL by varying the accretion rate (and other model parameters dependent on the accre- 
tion rate) in our specific spectral blazar model described in the previous section. In order to 
do so, we are starting out with parameters appropriate to reproduce the broadband spectrum 
of a typical FSRQ. The FSRQ phase constitutes an early phase of blazar evolution in which 
the central regions of the galaxy are rich in gas and dust, leading to a high accretion rate 
onto the central, supermassive black hole. At the same time, the circumnuclear material ef- 
ficiently reprocesses and scatters the accretion-disk radiation, leading to the observed strong 
optical emission lines in the broad line region and to a high energy density of the external 
soft photon field in the jet. 

Due to the limited supply of matter in the nuclear region, the average density of the 
circumnuclear material will gradually decrease, leading to a decreasing accretion rate and a 
decreasing reprocessing efficiency. We parametrize this evolutionary transition by a gradual 
decline of the accretion rate and assume that the reprocessing efficiency Tj-epr is proportional 
to the accretion rate. This should be a reasonable approximation for the long-term evolution 
on cosmological time scales, although the actual BLR density and the instantaneous accretion 
rate may be less strictly correlated on shorter time scales. Furthermore, we assume that the 
fraction of accreted energy channeled into the relativistic outfiow, fj, remains constant, as 
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suggested by the observed Lj vs. Lblr and Lj vs. Lnlr correlations (Rawling & Saunders 
1991; Celotti, Padovani & Ghisellini 1997; Xu, Livio & Baum 1999), and that the transverse 
extent of the jet, Rb, does not change dramatically during this evolution. The magnetic 
field is chosen to be a constant fraction of the equipartition magnetic field so that Ub = 
0.1 Me, which automatically leads to the required reduction of the magnetic field along the 
evolutionary sequence. 

The resulting sequence of broadband spectra, shown in Fig. 1, provides a quantitative 
explanation for the observed trend of luminosities and peak photon energies in the FSRQ — 
LBL HBL sequence (Fossati ct al. 1998). It docs not account for the spectral properties 
of specific blazars, which depend on variations of parameter values in different sources and, 
importantly, the angle between the direction of the jet axis and the direction to the observer. 
We note that the apparent minimum at ~ 10^° Hz in the FSRQ spectrum (corresponding to 
the Trepr = 0.1 curvc in Fig. 1) might be filled in by an additional contribution from Compton- 
scattered accretion disk radiation that enters the jet directly, which has been neglected in 
our semi-analytic model. 

Our estimate (Eq. [14]) for the average energy e* 2 x 10~^ of the dominant external 
soft photon field has interesting implications. Although one has to keep in mind that the 
analytical estimate for e* may be rather uncertain due to the uncertainty in the obscrvables 
pertaining to the SSC component, we found that significantly higher values of did not allow 
an acceptable representation of a typical FSRQ spectrum with our model. This probably 
indicates that the external soft photon field is dominated by optical line emission from the 
partially ionized broad-line regions. The external photon field peaking at is unlikely to be 
due to Thomson scattering off free electrons in a highly ionized circumnuclear environment, 
since in that case should reflect the peak energy of the optically thick accretion disk, which 
would be inconsistent with the evidence for thermal excess emission, peaking at ~ 2 x 10^^ Hz, 
observed in the low 7-ray state of 3C 279 by Plan et al. (1999) and with the peak frequency 
of ~ 3 X 10^^ Hz of the big blue bump seen in 3C 273 (Lichti et al. 1995) which is generally 
attributed to the emission from the optically thick accretion disk. 



4. Discussion and Summary 

The proposed evolutionary scenario provides a simple physical connection between dif- 
ferent blazar subclasses in terms of the depletion of the accretion flow onto supermassive 
black holes. After an early phase of rapid black-hole growth, the time scale for this depletion 
is given by the matter free-fall time tfj ~ 7 x 10^ mg T4 yr at the Bondi-Hoyle accretion 
radius, where mg is the mass of the central black hole in units of 10^ Mq and T = 10^^X4 K 
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is the average temperature of the circumnuclear material. This is generally shorter than 
the Hubble time, and gives a plausible time scale to explain a finite duration of blazar-like 
activity in the nuclei of elliptical galaxies. The scenario proposed here is also in accord with 
the X-ray faintness of giant elliptical galaxies known to contain supermassive black holes, 
which may be remnants of past blazar activity. The central sources in these galaxies have 
been successfully modeled with advection-dominated accretion flow solutions at very low 
accretion rates, I <g 1 (Fabian & Rees 1995; Di Matteo & Fabian 1997; Di Matteo et al. 
2000), suggesting a central region depleted after an earlier period of more efficient accretion. 

If this picture is correct, then there should also be objects in an earlier phase of rapid 
black- hole growth that are accreting at rates / ~ 1. Following the scenario outlined by 
Sanders et al. (1988), this earlier stage of blazar evolution would comprise merging galaxies, 
infrared luminous galaxies, and radio-quiet QSOs. Accretion flows near the Eddington limit 
might produce optically-thick accretion disks extending all the way to the innermost stable 
orbit. If the phenomenology of Galactic black-hole candidates can be scaled to supermassive 
black holes accreting at comparable values of I, then radio jet formation should be quenched 
in such a mode of accretion as is observed, e.g., in GRS 1915+105 (Mirabel et al. 1998; 
Feroci et al. 1999; BeUoni, Migliari, & Fender 2000), GX 339-4 (Fender et al. 1999; Corbel 
et al. 2000), and XTE J 1550-564 (Corbel et al. 2001). This idea is supported by recent 
theoretical work that hnks the existence and energetics of a central advection-dominated 
inflow-outflow system to the formation of relativistic jets (Blandford & Begelman 1999; 
Becker, Subramanian, & Kazanas 2001). The ASCA detections of relativistically broadened 
Fe Ka emission in MCG-6-30-15 (Tanaka et al. 1995) and IRAS 18325-5926 (Iwasawa et 
al. 1996) and the rapid X-ray variability seen in IRAS 13224-3809 and PHL 1092 (Boiler 
2000), suggest that at least in a number of cases an optically thick accretion disk docs extend 
towards the innermost stable orbit. In most cases the evidence is yet inconclusive, but a 
similar mode of accretion can generally not be ruled out. 

The formation of radio jets in blazars and radio galaxies, if related to the formation 
of an advection-dominated accretion mode in the inner portions of the accretion flow, is 
triggered by a decreasing Eddington ratio. The decline of / might be due to a combination of 
a decreasing accretion rate and an increasing black-hole mass. An implication of this scenario 
is that the masses of the central black holes in galaxies which host BL Lacs should, on average, 
be greater than the masses of black holes in galaxies which host FSRQs or optical QSOs. 
Moreover, subclasses at earlier stages in the blazar sequence should exhibit increasingly 
stronger cosmological evolution. Bade et al. (1998) flnd evidence for negative cosmological 
evolution of X-ray selected BL Lac objects, and Stickel et al. (1991) find evidence for positive 
cosmological evolution in a radio-selected BL Lac sample, in accord with this picture. Our 
work provides strong support for the suggestion by D'Elia & Cavaliere (2000) and Cavaliere 
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fe D'Elia (2001) along parallel lines, but on the basis of general energy considerations and a 
generic cosmological evolution scenario for quasars evolving into BL Lac objects. We have 
combined the basic idea of such a cosmological evolution with a self-consistent spectral blazar 
model and demonstrated that the spectral properties of the different blazar subclasses can 
be consistently explained by such an evolutionary sequence. The blazar model adopted here 
is similar to the one previously used by Ghisellini et al. (1998) — but with the important 
addition of particle escape and adiabatic losses — who had argued that the different blazar 
subclasses could be unified by an increasing level of external radiation to the soft radiation 
field in the emitting region. 

To summarize, we have proposed an evolutionary scenario linking FSRQs, LBLs, and 
HBLs through gradual depletion of the circumnuclear environment of a supermassive black 
hole. As the accretion power declines with time, less gas and dust is available to reprocess 
accretion-disk radiation and produce an external Compton-scattered component in the blazar 
jet. Analytic formulae are given to estimate the relevant model parameters for FSRQs. We 
have calculated a spectral sequence for different blazar subclasses that are related through 
a one-parameter model defined by the optical depth of the circumnuclear matter. This 
scenario links radio-loud AGNs with progenitor merger galaxies and QSOs that constitute 
the high Eddington-ratio limit of the evolutionary sequence. In analogy with the observed 
anti-correlation between radio and soft X-ray activity in some Galactic black-hole candidates, 
these progenitor sources would accrete with large Eddington ratios, resulting in an optically- 
thick accretion disk that quenches relativistic jet production. 

We thank the anonymous referee for a very careful review and helpful and constructive 
comments. The work of MB is supported by NASA through Chandra Postdoctoral Fellow- 
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Fig. 1. — One-parameter model sequence of broadband blazar spectra that reproduces the 
trend in the spectral energy distributions of FSRQs, LBLs, and HBLs. Starting with typical 
FSRQ parameters, 71 = 100, 72 = 10^ s = 2.2, D = 28, Rb = 3 x 10^^ cm, B = 2.4 G, 
= lO"'^ ergs s^^, Trcpr = 0.1, = 10^^, fj^t = 10^'^, Rsc = 1-2 x 10"*^^ cm, the sequence is 
generated by reducing Tj-cpr and assuming that the accretion disk luminosity is proportional 
to Trepr- The comoviug magnetic field B = 0.3 Bgp, where Bgp is the equipartition magnetic 
field with electrons. 



